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(57) Abstract: The invention relates to a method for determining the position of a sensor element (300), according to which a 
magnetic alternating field (210) emitted by al least one field generating unit (200) is measured. The position of the sensor element 
(300) is determined on the basis of a signal received in the sensor element (300). The inventive method is further characterized in that 
interference fields (410) are calculated, preferably to a first approximation, said interference fields being caused by eddy currents 
(420) produced in electrically conductive objects (400). The position that can be determined on the basis of the signal received in 
the sensor element (300) is corrected on the basis of the calculated interference fields (410). 
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(57) Zusammenfassung: Die Erllndung betrifft zunachst ein Verfahren zur Bestimmung der Position eines Sensorelememes (300), 
mil Hilfe dessen ein von mindestens einer Feldgeneratoreinheit (200) ausgesendetes magnetisches Wechselfeld (210) gemessen 
wird, wobei aufgrund eines im Sensorelemeni (300) empfangenen Signals die Position des Sensorelememes (300) bestimmt wird. 
Erfindungsgemass ist vorgesehen, dass, vorzugsweise in erster Naherung, Storfelder (410) berechnet werden, welche aufgrund von 
in elektrisch leiienden Objekien (400) erzeuglen Wirbelstromen (420) entstehen. Des weiteren ist vorgesehen, dass die Position, 
welche ausgehend vom im Sensorelemeni (300) empfangenen Signal bestimmbar ist, aufgrund der berechneten Storfelder (410) 
konigiert wird. 
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Verfahren zur Bestimmung der Position eines Sensorelementes 

Die vorliegende Erf indung betrifft ein Verfahren nach dem 
5 Oberbegriff des Patentanspruchs 1, eine Anwendung des 
Verfahrens, eine Vorrichtung zur Durchfuhrung des 
Verfahrens sowie ein Computerprogrammprodukt . 

Bei zahlreichen technischen und medizinischen Verfahren ist 
10 die genaue Kenntnis der Position eines bestimmten 

Gegenstandes von entscheidender Bedeutung. Wahrenddem in 
der Medizin die Position von einzelnen Gewebeteilen - 
beispielsweise eines Tumors, der zur Zerstorung oder zur 
Wachstumsbegrenzung bestrahlt werden soil - bestimmt werden 
15 muss, ist die Positionserf assung zur Eingabe in ein 
Computersystem, beispielsweise fur „Cyber Space*- 
Anwendungen, von allgemeiner Bedeutung. Eine solche 
Positionserf assungs- bzw. Positionseingabeeinheit wird in 
diesen Anwendungen auch etwa als dreidimensionale Maus 
20 bezeichnet. 

Eine bekannte Vorrichtung bzw. ein bekanntes Verfahren zur 
Bestimmung der Position ist in der internationalen 
Patentanmeldung mit der Verof f entlichungsnummer WO 97/36192 

25 der gleichen Anmelderin beschrieben. Gemass der bekannten 
Lehre ist es vorgesehen, ein Wechselfeld mit Hilfe einer 
Feldgeneratoreinheit aufzubauen, wobei je nach der Anzahl 
Freiheitsgrade eines Sensorelementes, dessen Position 
bestimmt werden soli, mehrere Wechself elder einander 

30 uberlagert werden. Mit Hilfe einer Verarbeitungs- und 

Steuereinheit , welche die Feldgeneratoreinheit einerseits 
steuert, anderseits die vom Sensorelement empfangenen 
Signale verarbeitet, wird die Position und gegebenenf alls 
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die Lage der Sensoreinheit bestimmt. Der Inhalt der 
vorstehend genannten Verof f entlichung bildet diesbezuglich 
integrierender Bestandteil dieser Beschreibung. 

5 Es hat sich gezeigt, dass bei einer magnetfeldbasierten 
Ortung, wie sie beispielsweise bei der bekannten Lehre 
gemass WO 97/36192 angewandt wird, in benachbarten, 
elektrisch leitenden Objekten Wirbelstrdme erzeugt werden. 
Diese fuhren zu Verzerrungen des ursprunglichen 
10 magnetischen Wechselfeldes und damit zu systematischen 

Fehlern. Das bedeutet, dass, wenn Position und Ausrichtung 
von Sensorelementen im verzerrten Wechselfeld so bestimmt 
werden, als ware kein elektrisch leitendes Objekt zugegen, 
die gewonnenen Werte systematise!! verfalscht sind. 

15 

Ein Verf ahren zur Kompensation von Storef f ekten, welche 
durch leitende Objekte hervorgeruf en werden, ist unter dem 
Namen „ Distortion Mapping" bekannt. Dieses Verf ahren wird 
beispielsweise im Aufsatz mit dem Titel „ Calibration of 

20 Tracking Systems in a Surgical Environment* (Birkfellner 

et. al., IEEE Trans Med Imaging, Vol- 17(5), Seiten 737 bis 
742, 1998) beschrieben. Beim bekannten Verf ahren wird die 
Position und die Orient ierung eines Sensorelementes 
ebenfalls mit Hilfe eines Positionsmesssystems, das auf 

25 einer magnetfeldbasierten Ortung beruht, vorgenommen, wobei 
zur Kompensation von Storef fekten ein zweites 
Positionsmesssystems vorgesehen ist, welches durch 
elektrisch leitende Objekte nicht beeinf lussbar ist. Die 
Differenz zwischen den mit den beiden Positionsmesssystemen 

30 ermittelten Posit ionen und Orient ierungen werden in der 

Folge zur Korrektur der mit Hilfe des magnetfeldbasierten 
Positionsmesssystems bestimmten Position und Orientierung 
verwendet. Das bekannte Verf ahren weist jedoch den Nachteil 
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auf, dass, zur Erreichung einer hohen Genauigkeit, die 
Positions- und Orient ierungsdif f erenz an moglichst vielen 
Punkten gemessen werden muss- Urn weitere Punkte zu 
erhalten, muss zusatzlich ein aufwendiges 
5 mterpolationsverfahren angewendet werden. Der sehr hohe 
Auf wand wird insbesondere durch das folgende Beispiel 
deutlich: Soil ein Volumen von 1 m 3 vermessen werden, wobei 
dies in den drei Achsen alle 10cm und in zehn verschiedenen 
Orientierungswinkeln geschehen soli, so erhalt man 10' 000 
10 Punkte. Daruber hinaus ist das erwahnte zweite 
Positionsmesssystem erforderlich. 

Des weiteren ist ein Verfahren zur Kompensation von 
Storeffekten bekannt, bei dem Magnetf elder durch gepulste 

15 Gleichstromf elder erzeugt werden, wobei die Kompensation 
von Wirbelstromef fekten dadurch vorgenommen werden, indem 
Magnetfeldmessungen erst nach dem Abklingen der im 
Messsignal enthaltenen Wirbelstromanteile durchgefiihrt 
werden. Weiterffthrende Erlauterungen zum bekannten 

20 Verfahren konnen den Druckschrif ten US-5 453 686 und US-5 
767 669 entnommen werden. Es hat sich gezeigt, dass die 
Genauigkeit der ermittelten Resultate ungenugend ist. 
Insbesondere ist die Kompensation unvollstandig, wenn die 
Zerfallszeiten der Wirbelstromanteile die Pulszeit zwischen 

25 zwei Gleichstrompulsen ubersteigt. Zwar kann dem durch 

Verlangerung der Pulszeit begegnet werden, doch fuhrt dies 
zu einer unerwunschten geringeren Messrate. Des weiteren 
lasst sich die bekannte Kompensationsmethode nicht bei auf 
magnetischer Ortung basierenden Positionsmesssystemen 

30 anwenden, die magnetische Wechselfelder erzeugen. 

Der vorliegenden Erf indung liegt daher die Aufgabe 
zugrunde, ein Verfeihren anzugeben, das eine verbesserte 
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Bestimmung der Position und/oder der Lage eines 



Sensorelementes ermoglicht. 

Diese Aufgabe wird durch die im kennzeichnenden Teil des 
5 Patentanspruchs 1 angegebenen Massnahraen gelost. 

Vorteilhafte Ausgestaltungen der Erfindung, eine Anwendung 
des Verf ahrens, eine Vorrichtung zur Durchfiihrung des 
Verf ahrens sowie ein Computerprogrammprodukt sind in 
weiteren Anspruchen angegeben. 

10 

Mit dem erf indungsgemassen Verfahren wird es ermoglicht, 
den Einfluss von leitenden Objekten zu eliminieren, 
zumindest aber erheblich zu reduzieren. Des weiteren ist 
dieses Verfahren allgemeiner und genauer als die bekannten 
15 verfahren. Schliesslich kann der geometrieabhangige Teil 
der Berechnungen im Sinne einer Systemkalibrierung vor dem 
eigentliche Einsatz des Positionsmesssystems vorgenommen 
werden. 

20 Die Erfindung wird nachfolgend anhand von Zeichnungen bei- 
spielsweise naher erlautert, Dabei zeigen 

Fig- 1 eine bekannte Anordnung, bestehend aus 



25 



Feldgeneratoreinheit, Sensorelement xxnd 
Verarbeitungs- und Steuereinheit , in 
schematischer Darstellung, mit einem elektrisch 
leitenden Objekt, 



Fig. 2 



ein elektrisch leitendes Objekt und 



30 



Fig. 3 



ein Flussdiagramm mit einigen Verf ahrensschritten 
des erf indungsgemassen Verf ahrens. 
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In Fig- 1 ist eine bekannte Anordnung, bestehend aus einer 
Feldgeneratoreinheit 200, einem Sensor element 300 und einer 
Verarbeitungs- und Steuereinheit 100, dargestellt. Die 
Verarbeitungs- und Steuereinheit 100 ist jeweils uber 
5 Leitungen mit der Feldgeneratoreinheit 200 einerseits und 
dem Sensorelement 300 anderseits verbunden. Wahrend sich 
die Feldgeneratoreinheit 200 vorzugsweise an einem 
bekannten Ort befindet - was bedeutet, dass die Koordinaten 
x, y, z, inkl. der Ausrichtung im Koordinatensystem, 

10 bekannt sind - kann das Sensorelement 300 beliebig bewegt 
werden bzw. eine beliebige Position und Orient ierung 
einnehmen. Es wird darauf hingewiesen, dass es denkbar ist, 
wie bereits aus der WO 97/36192 bekannt ist, dass das 
Sensorelement 300 ort sf est und die Feldgeneratoreinheit 200 

15 frei, d.h. im Rahmen der zur Verfugung gestellten 

Verbindungsleitung zur Verarbeitungs- und Steuereinheit 
100, beweglich ist. Des weiteren ist auch ohne weiteres 
denkbar, dass die Verarbeitungs- und Steuereinheit 100 in 
mehreren Funktionseinheiten realisiert ist, wie zum 

20 Beispiel, dass die Steuereinheit zur Steuerung der 

Feldgeneratoreinheit 200 in einem Funktionsblock und die 
Verarbeitungseinheit, in der die eigentliche 
Positionsberechnung der Position des Sensorelementes 300 
vorgenommen wird, in einem anderen Funktionsblock 

25 realisiert ist. Diese Abwandlungen von der in Fig. 1 
dargestellten Anordnung haben keinen Einfluss auf die 
Anwendbarkeit des erf indungsgemassen Verfahrens. Das 
Gleiche gilt auch fiir Ausfuhrungsf ormen, bei denen mehrere 
Feldgeneratoren an unterschiedlichen Stellen vorgesehen 

30 sind, wie dies beispielsweise bei der Lehre gemass WO 
97/36192 der Fall ist. 

Mit 400 ist in schematischer Darstellung ein elektrisch 
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leitendes Objekt stellvertretend fur diejenigen Gegenstande 
dargestellt, welche die magnetische Ortung des 
Sensorelementes 300 stort, indem im Objekt 400 Wirbelstrome 
420 erzeugt werden, aufgrund derer ein dem Wechselfeld 210 
5 fiberlagerten Storfeld 410 entsteht. 

Bevor auf das erf indungsgemasse Verfahren weiter 
beschrieben wird, werden im folgenden zunachst allgemeine 
Zusamraenhaiige bzw. Vorgehensweisen bei der 
10 magnetfeldbasierten Ortung erlautert. 

Wie bereits erwahnt wurde, wird bei der magnetfeldbasierten 
Ortung, welche auch etwa als magnetische Ortung bezeichnet 
wird, die Position und/oder Ausrichtung von einem oder 

15 mehreren Sensorelementen 300 relativ zu einem oder mehreren 
Feldgeneratoreinheiten 200 bestimmt. Die Position r Si und 
Orientierung n Si der Sensorelemente S, kann durch Losung des 
folgenden Gleichungssystems bestimmt werden, wenn man 
voraussetzt, dass Position r GJ und Orientierung H CJ der 

20 Feldgeneratoreinheiten G } bekannt sind: 

F & = F{r St9 n Si ,r Gj ,n Gj ) (1) 

Mit i ist das i-te Sensorelement und mit j die j-te 
25 Feldgeneratoreinheit gemeint. F ist dabei eine vom 
Magnetfeld abhangige Messfunktion von meistens einer 
Komponente des Magnetfelds B{x,y,z,t) (z.B. die induzierte 
Spannung in einer Sensorspule) . F kann naturlich auch eine 
Funktion von vielen zusammengebauten Sensoren in einem 
30 Sensorelement sein, welches mehrere oder alle Komponenten 
gleichzeitig misst. 
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Nach Art der Losung dieses Gleichungs systems konnen 
magnetische Positionssysteme in zwei Klassen unterteilt 
we r den: 

5 

I. Das Gleichungssystem wird invert iert, d.h. die 
Sensorelementpositionen konnen aus den gemessenen 
Magnet feldern berechnet werden: 

10 r Sl ^f r (F 9 )undn Si =/ B (^) (2) 

Da die Inversion des Gleichungssystems nur in sehr 
speziellen Fallen moglich ist, kann man durch Naherung 
versuchen, die Feldgleichungen in eine invert ierbare Form 
15 zu bringen. 

II. Das Gleichungssystem wird durch Optimierung gelost, 
d.h. die Sensorelementpositionen werden so lange variiert, 
bis die nach Gleichung 1 berechneten Werte F 9 am besten 

20 mit den gemessenen Werten F g M ubereinstimmen. Eine 

mogliche Methode ware ein Chi 2 -Fit nach Levenberg- 
Marquardt. Dabei werden die Sensorpositionen r Si und n Si so 

lcinge variiert bis 



25 



J (AFg ) 



minimal wird. Fur weiterfuhrende Angaben zur Methode nach 
Levenberg-Marquardt wird stellvertretend auf die 
Druckschrift mit dem Titel „Numerical Recipies in C w (W. H. 
30 Press, S. A. Teukolsky, W. T. Vetterling und B. P. 

Flannery; Cambridge university Press; 1994) verwiesen. 
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Es ist auch eine Kombination beider Losungsansatze moglich. 

Da die vom Sensorelement Si gemessene Grosse F 9 nur von der 

5 relativen Lage vom Sensorelement S ± und der 

Feldgeneratoreinheit abhangt, sind die Rollen von 
Sensorelement und Feldgeneratoreinheit in alien 
magnetischen Positionsmesssystemen austauschbar . 

10 Werden zeitlich variierende magnetische Felder bentttzt, so 
erzeugen diese - wie bereits erwahnt wurde - in 
benachbarten, elektrisch leitenden Objekten 400 
Wirbelstrome 420. Diese fuhren zu Verzerrungen des 
ursprunglichen magnetischen Wechself eldes 210 und damit zu 

15 systematischen Fehlern bei der Positionsbestimmung. Das 
bedeutet, dass, wenn Position \ind Ausrichtung von 
Sensorelementen im verzerrten Wechselfeld so bestimmt 
werden, als ware kein elektrisch leitendes Objekt 400 
zugegen, die gewonnenen Werte systematisch verfalscht sind. 

20 

Um das Verfahren der magnetischen Ortung auch in der 
Nachbarschaft von elektrisch leitenden Objekten 400 frei 
von durch diese hervorgeruf enen Fehlern anwenden zu konnen, 
werden erf indungsgemass die Wechself eldverzerrungen und ihr 
25 Einfluss auf die Bestimmung von Sensorelementposition und 
Sensorelementausrichtung bestimmt. Damit konnen die 
auftretenden systematischen Fehler korrigiert werden, womit 
die Genauigkeit der Position und/oder der Ausrichtung 
erheblich verbessert werden kann. 

30 

Bei einer vorgegebenen Messanordnung hat man feststellen 
konnen, dass ein Mess fehler von 4cm mit dem 
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erfindungsgemassen Verfahren auf weniger als 1,5mm 
reduziert werden konnte. 

Grundsatzlich lassen sich diese Korrekturen auch mit Hilfe 
5 der Technik der f initen Elemente und den Gleichungen der 
Elektrodynamik finden. Die bevorzugte Ausfuhrungsf orm des 
erfindungsgemassen Verfahrens zeichnet sich gegenuber der 
Methode der f initen Elemente weiter dadurch aus, dass eine 
massive Reduktion der Positionsberechung erreicht wurde, da 
10 vieles im Sinne einer Systemkalibration im voraus berechnet 
werden kann. 

Im folgenden wird das erf indungsgemasse Verfahren 
erlautert, wobei der Einfachheit halber zunachst davon 

15 ausgegangen wird, dass das Objekt 400 aus einer elektrisch 
leitenden Platte, d.h. aus einer flachen und begrenzten 
Flache, besteht. Objekte 400, die in der Richtung einer 
gedachten Linie von der Feldgeneratoreinheit 200 zum Objekt 
400 eine relevante Ausdehnung (Tiefe) besitzen, konnen 

20 ebenfalls mit dem erfindungsgemassen Verfahren behandelt 
werden. Fur diesen Zweck wird die der Feldgeneratoreinheit 
200 zugewandte Seite durch eine Vielf lachen-Struktur 
angenahert. Dies ist zulassig, da die Wirbelstrome 420 nur 
geringfugig in die Oberflache eindringen. Die Tiefe eines 

25 dreidimensionalen Objektes 400 ist deshalb nicht relevant. 
Bei der mathematischen Bestimmung des Storfeldes 410 wird 
daher das Objekt 400 durch eine Vielf lachenstruktur im 
vorstehend genannten Sinne angenahert. 

30 Im folgenden werden Uberlegungen betreffend die Berechnung 
von Feldverzerrungen bei einer leitende Platte dargestellt. 
Die daraus resultierenden Ergebnisse konnen analog bei 
allgemeineren Objektformen angestellt werden. 
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Bef indet sich in einem zeitlich verander lichen Magnetfeld 
B 0 (x 9 y 9 z 9 t) ein elektrisch leitendes Objekt 400 , so werden 
Wirbelstrome 410 (Fig. 1) in die Oberflache des Objektes 
5 400 induziert. Diese Wirbelstrome 410 verursachen ein 

weiteres Magnetfeld B\x 9 y 9 z 9 t) , welches dem ursprunglichen 
Feld B 0 (x 9 y 9 z 9 t) uberlagert ist und zu einem resultierenden 
Feld B^ s (x 9 y 9 z 9 t) fuhrt. B^ s (x 9 y 9 z 9 t) ist gegenQber dem Feld 
B 0 (x 9 y 9 z 9 t) verzerrt. Dm dieses verzerrte Feld berechnen zu 
10 konnen, ist es notwendig, das induzierte Wechselfeld 
B\x 9 y 9 z 9 t) zu kennen, Es lasst sich aus dem Biot-Savart 
Gesetz der Elektrodynamik (Gleichung 4) ein Feld B x (x 9 y 9 z 9 i) 

berechnen, welches gut genug B\x 9 y 9 z 9 t) beschreibt, falls 
der ortliche und zeitliche Verlauf der Wirbelstrome 410 im 
15 Objekt 400 in N verschiedenen punkt formigen Stromelementen 
bekannt ist: 

20 wobei P R ={x 9 y 9 z) ein Punkt im Raum ist und der Vektor r vom 
Stromelement zum Punkt P R zeigt. Vorfaktoren konnen, falls 
notwendig, eingeffihrt werden und/oder konnen im Bet rag des 
Faktors As(t) enthalten sein. Eine detailliertere Berechnung 

von B x (x 9 y 9 z 9 t) durch Einfuhrung von Langsstromen oder 
25 Flachenstromen etc. anstatt Punktstromen ist moglich. 

Dadurch wurde die Schreibweise von Gleichung 4 allerdings 
geandert werden. In den meisten Fallen kann die Gleichung 4 
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hingegen so \ibernommen werden, wie sie oben aufgefuhrt ist # 
falls N gross genug gewahlt wird. 

Die Feldverzerrungen werden also in zwei Schritten 
5 berechnet. Der erste Schritt ist die Bestimmung der 

Wirbelstrome 410 und der zweite Schritt ist die Berechnung 
eines durch die Wirbelstrome 410 erzeugten Storfelds 
B x (x 9 y 9 z 9 t), welches gut genug das Storfeld B 9 (x 9 y 9 z 9 t) 
beschreibt. 

10 

Fig- 2 zeigt das Objekt 400 , das fur die Bestimmung des 
Storfeldes 410 in eine Vielf lachenstruktur, bestehend aus 
beliebig vielen Segmenten, aufgeteilt ist. Ausgehend von 
dieser Auf teilung und diversen weiteren Annahmen werden 
15 zunachst die Wirbelstrome berechnet. 

Die Wirbelstrome fliessen an der Oberflache des Objektes 
400 mit einer fur die Theorie "unwesentlichen" 
Eindringtiefe. Urn das oben erwahnte Storfeld B'(x 9 y 9 z 9 t) 

20 hinreichend genau zu berechnen, reicht es aus, den 
zeitlichen Stromverlauf in einigen Punkten auf der 
Oberflache des Objektes 400 zu kennen. Die Anzahl der 
Punkte hangt von der verlangten Genauigkeit ab. Die 
Wirbelstrome werden also in Punkten berechnet, welche auf 

25 Oder in der Nahe der Oberflache des Objektes liegen. 

In einem ersten Schritt wird das Objekt in N beliebig 
geformte Segmente unterteilt, die sinnvollerweise (aber 
nicht notwendigerweise) das ganze Objekt abdecken. Die 
30 Segmente werden nachfolgend zur eindeutigen Unterscheidung 
mit S;{0£i<N-l} bezeichnet, wobei i als Index verwendet 
wird- 
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In einem zweiten Schritt wird pro Segment ein Stutzpunkt /> 
gewahlt. Es ist sinnvoll aber nicht zwingend, gleich viele 
Stutzpunkte wie Segmente zu definieren und diese eindeutig 
5 zu den Segmenten zu zuordnen. Im folgenden wird der 

Einfachheit halber davon ausgegangen, dass N Segmente S t 
mit je einem eindeutig zugeordneten Stutzpunkt P t definiert 
werden. Die Stromdichte 7,(0 im Stutzpunkt P t eines jeden 
Segments S, wird mittels der folgenden Formel berechnet: 

10 

wobei t!j(t) die Stromdichte des Wirbelstroms 7^(0 ist, 

welche durch die Flussanderung des Feldes von B 0 (x,y,z,t) im 
15 Segment S f verursacht wird und durch den Stutzpunkt P t Oder 
in der Umgebung des Stutzpunktes P t fliesst. Die Berechnung 
der einzelnen Wirbelstrome I g (t) ist im nachsten Abschnitt 
beschrieben. Zunachst gilt: 



20 




(5a) 



mit X als Richtungsvektor (oder eine dazu fast kolineare 
Richtung) der Stromlinie durch den Stutzpunkt P f im 
Stutzpunkt P s und mit A s die Querschnittsf lache der 
25 Stromline ist, wobei: 

A,=xr 2 (h) (5b) 



BNSOOCIO: <WO 0208793A1J_> 



WO 02/08793 



PCT/CH0 1/00431 



- 13 - 

mit r = Radius der kreisrunden Querschnittsf lache; 
h = Eindringtiefe; 

Sind nun die Stromdichten ^(/) bekannt, so ist auch 
5 B } (x 9 y,z 9 t) als Storfeld verursacht durch B 0 (x 9 y 9 z 9 t) 

berechenbar - dazu kann %{tyA(S,) direkt in die Gleichung 4 
eingesetzt werden, wobei A(S t ) die Flache des Segments S A 

ist. B x {x 9 y 9 z 9 t) kann in den meisten Fallen als B'(x,y 9 z,t) 
betrachtet werden. Als Effekt zweiter Ordnung kann an 

10 dieser Stelle B x (x 9 y 9 z 9 t) als Urfeld eingesetzt werden, urn 

wiederum Wirbelstrome fur ein zweites Storfeld B 2 (x 9 y 9 z 9 t) zu 
berechnen (Einf luss der Wirbelstrome auf einander) , welches 
mit B 0 (x 9 y 9 z 9 t) und B x (x 9 y 9 z 9 t) uberlagert ist. In zweiter 
Naherung ware also B'(x 9 y 9 z 9 t) gleich der Summe aus B x {x 9 y 9 z 9 t) 

15 und B 2 (x 9 y 9 z 9 t) - diese iterative Vorgehensweise kann ffir 
Effekte beliebiger Ordnung fortgesetzt werden. Es ist 
jedoch in den meisten Anwendungen der Effekt erster Ordnung 
hinreichend genau. 

20 Ein einzelner Wirbelstrom I 9 (t) ist eine Stromlinie, die 
durch den Stutzpunkt P ( fliesst und durch die zeitliche 
Flussanderung des Feldes durch den Stutzpunkt P. verursacht 
wird. Um I y (t) zu berechnen sind die Induktivitat L 9 , sein 
Ohm'scher Widerstand R v und die zeitliche Flussanderung 

25 notwendig. Sind diese Grossen bekannt, so ist I g (t) 

dt 

durch die Losung der bif f erentialgleichung 
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gegeben. In vielen Fallen mag B Q (x,y,z 9 t) zeitlich periodisch 
sein oder sogar harmonisch oszillieren, dies ist aber fur 
5 die Gultigkeit der erf indungsgemassen Methode keine 
Notwendigkeit . 

Die Induktivitat L v und der Ohm'sche Widerstand R 0 sind 
durch die geometrische Form des Wirbelstroms I g (t) gegeben 

d<!> 

10 und die Flussanderung — L durch das Feld B 0 (x,y,z 9 t) an der 

at 

Stelle Pj zusammen mit der Flache des Segments 5 y . In den 
folgenden Schritten wird zuerst die Form des Wirbelstroms 
beschrieben, und anschliessend daraus die Induktivitat L i} 
und der Ohm'sche Widerstand R g berechnet. 



15 



20 



Man stelle sich vor, dass eine einzelne Magnetf eldlinie B 
durch eine kleine Flache dA urn einen Punkt P auf dem 
Objekt dringt. In diesem Fall waren die induzierten 
Stromlinien in der Nahe vom Punkt P kreisformig und am 
Rande des Objektes 400 wurden sie der Begrenzung folgen, 
also die Form der Umrandung des Objekts haben. Die Form 
eines beliebigen Wirbelstroms ist ein Strom entlang einer 
Hohenlinie einer Flache, welche die Potentialgleichung 



25 A^ = |j^ = 0 (7) 

2 * Y dx 1 dy 1 



erffillt, wobei <p(x 9 y) das Potential darstellt. Die 
Randbeding\ingen zur eindeutigen Losung von Gleichxing 7 
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(ermitteln von <p(x,y)) sind aus Figur 2 zu entnehmen 
(namlich <J> 0 am Rand des Objektes und ^ im Punkt P j# <|> 0 <> 
ty x ) . Diese Potentialgleichung ist am besten numerisch zu 
losen. Mit der Form des Wirbelstroms und der Eindringtiefe 
5 h des Stroms in das Material ist die einzelne Stromlinie 
als Leiterschleife mit einem ringformigen 
Materialquerschnitt mit dem Durchmesser der Eindringtiefe 
betrachtbar (andere sinnvolle 

Materialquerschnitt sgeometrien sind natfcrlich denkbar, 
10 Sndern aber an den Berechnungen nichts wesentliches) . 

Der Ohm'sche Widerstand der Leiterschleife wird somit zu 



(8 ) 



15 



25 



is)" 



mit R ■ ohmscher Widerstand [ SI ] 

1 = Lange der Leiterschleife t m ] 
h a Eindringtiefe des Stroms [ m ] 

p a spezifischer elektrischer Widerstand des Materials 
20 [ fim ] 

Die Induktivitat der Leiterschleife ist gegeben durch 



2W 



und kann numerisch berechnet werden, wobei 



W = f -H-^dV (10) 
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die gespeicherte Energie im vom Leiter erzeugten Magnetfeld 
bedeutet, falls ein Strom i im Leiter fliesst. Es gibt 
"beliebig" viele andere Naherungsformeln, welche 
5 Gleichungen (8) , (9) und (10) ersetzen konnten und ahnliche 
Resultate liefern. 

Fur die Berechnung des Flusses ®j(t) aus dem Feld B Q {x 9 y,z y t) 
gilt : 

10 

Oj(t) = B 0 (x 9 y 9 z 9 t)^Aj (11) 

wobei B 0 (x 9 y 9 z y t) das ungestorte Feld an der Stelle P j ist 

und Aj die Flachennormale des Segments Sj ist, mit dem 
15 Betrag der Flache des betreffenden Segments. Formel 11 ist 
eine Naherungsformel fixr die allgemein gultige Formel 

Q>j(t) = lB Q (x f y,z y t)dA (lla) 

und darf angewendet werden, wenn B uber Aj hinreichend 
homogen ist (so z.B. fur kleine Flachen Aj) • An dieser 
20 Stelle soli erwahnt werden, dass die Flache A } des 

Segments S } nicht in alien Fallen vollstandig innerhalb der 
Leiterschleife I & {t) liegt - Korrekturen konnten 
diesbezuglich angebracht werden, sind aber in der Regel 
nicht notwendig. 

25 

Die teilweise intensiven Berechnung der Induktivitaten L y 
und der Ohm'schen Widerstande R & konnen mittels der 
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Formeln (8), (9) und (10) im voraus berechnet werden, da 
diese allein von der Geometrie und des Materials des 
Objektes abhangig sind - dies im Sinne einer 
Systemkalibration. Einsetzen und Losen der Formeln (4) und 

5 (5) , um 5,(jc,^,z,/) zu berechnen, konnen zu einem anderen 
Zeitpunkt erfolgen, insbesondere wenn das Erregerfeld 
B 0 (x 9 y 9 z 9 t) bekannt ist. Iterationen, wie z.B. Einsetzen des 
Feldes B x (x 9 y 9 z,t) zur Berechnung eines Feldes B 2 (x 9 y 9 z 9 t) 
usw., sind moglich. Solche Iterationen konnten aber auch im 
10 voraus erfolgen, indem man sie in die Induktivitaten L fJ 

einbindet. Dies macht aber erst einen Sinn, wenn 
regelmassig Korrekturien hoherer Ordnung benotigt werden. 

In der Praxis werden bei der magnetischen Ortung Lage und 

15 Ausrichtung eines oder mehrerer Sensorelemente 300 (Fig. 1) 
in einem Magnet f eld bestimmt, das von einem oder mehreren 
Feldgeneratoreinheiten 200 erzeugt wird. In dem verwendeten 
Koordinatensystem ist die Position der Feldgeneratoreinheit 
200 oder den Feldgeneratoreinheiten bekannt. Im Fall von 

20 magnetischen Wechself eldern erzeugen benachbarte elektrisch 
leitende Objekte 400 Feldverzerrungen durch in den Objekten 
400 induzierte Wirbelstrome 420. Das erf indungsgemasse 
Verfahren zur Korrektur dieser Verzerrungen, dessen 
theoretische Grundlagen oben angegeben wurden, wird 

25 folgendermassen angewendet: Die Lage der elektrisch 

leitenden Objekte 400 im oben erwahnten Koordinatensystem 
ist bekannt oder wird durch Vermessung bestimmt. Die 
Objektkoordinaten werden in einem Computerprogramm, das zur 
Berechnung der Wirbelstrome 420 und der daraus 

30 resultierenden Feldverzerrungen benutzt wird, derartig 
eingegeben, dass die in den oben angegebenen Formeln 
verwendeten Ortskoordinaten im durch die 
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Feldgeneratoreinheit 200 definierten Koordinatensystem 
definiert sind. Mit dem Computerprogramm wird dann das 
durch die Wirbelstrome 420 erzeugte Storfeld berechnet. Bei 
Berucksichtigung der Wirbelstrome 420 andert sich das 
5 Gleichungssystem 1 wie folgt : 

p 

iV=F(r 5/ ,w v V^ (12) 

wobei F 0 ' die durch den Wirbelstrome 420 des Objektes k 

10 erzeugte Stoning darstellt. P ist die Anzahl der Objekte. 
Wie diese Korrektur angewendet wird, hangt von der Art des 
magnetischen Positionsmesssystems ab. 

I. Bei Systemen, die auf Gleichung 2 beruhen, werden die 
15 Messwerte iterativ korrigiert, d.h. man berechnet zuerst 

die ungestorte Losung nach Gleichung 2. Mit der gefundenen 
Position des Sensorelementes 300 kann der Korrekturterm F 1 

berechnet und von den Messungen F 9 M abgezogen werden. Mit 

den korrigierten Messungen wird wieder eine Position 
20 berechnet. Dieser Algorithmus wird fortgesetzt, bis die 
Variation der berechneten Positionen unter gewissen 
Toleranzschwellen liegen. 

II. Bei Systemen, die auf Gleichung 3 beruhen, muss der 

25 Ldsungsalgorithmus nicht geandert werden. In der Chi 2 -Summe 
werden anstatt des Modells zur Berechnung der 
wirbelstromf reien Magnetfeldes F y nach Gleichung 1 das 
Magnetfeld mit Wirbelstromkorrekturen nach Gleichung 12 
benutzt . 
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III. Unter bestimmten Voraussetzungen kann es auch moglich 
sein, das Gleichungssystem 12 zu invertieren, was dann zu 
einer Losung entsprechend Gleichung 2 fuhrt. 

5 Fig. 3 zeigt, in vereinfachter Darstellung, ein 

Struktogramm eines nach dem erf indungsgemassen Verfahren 
arbeitenden Computerprogramms . Die einzelnen 
Verarbeitungsschritte wurden bereits anhand der Pig. 1 und 
2 ausfuhrlich erlautert. 

10 Das erfindungsgemasse Verfahren kann auch vorzuglich fur 

Objekte roit Offnungen (Locher) angewendet werden, wobei die 
Anzahl L der Offnungen beliebig sein kann. Dabei muss fur 
die schon beschriebene Losungsmethode zunachst die 
Randbedingung der Potent ialgleichung (7) an den Randern der 

15 offnungen gleich dem Potential <|> 0 am Rande des Objektes 
sein. Des weiteren kommen zu der beschriebenen Methode 
weitere N mal L (N=Anzahl Stutzpunkte und L=Anzahl 
Offnungen) Stromlinien I ik hinzu, die in der Summe (5) noch 
zu addieren sind (k wandert von 1 bis L) . 

20 Die zusatzlichen Wirbelstromlinien I ilt sind einzeln analog 
zu den Wirbelstromlinien berechenbar, d.h. Losen der 
Potent ialgleichung (7) fur die Form des Stroms und 
berechnung der Induktivitat und des Widerstands nach 
Gleichungen (8) und (9) . Bei der Potentialgleichung (7) ist 

25 allerdings zu beachten, dass die Randbedingung nicht im 
Punkt Pj n lautet, sondern am Rand der Offnung k" . Bei 
grossen Offnungen ist eventuell Formel 11a fur die 
Berechnung des Flusses zu verwenden anstatt die 
Naherungsf ormel 11. 

3 0 Einzelne Leiterschleif en konnen auch mit dieser Methode 
berechnet werden, da die vorher erwahnten Offnungen 
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beliebig nahe an der Umrandung der zu berechnenden Objekte 
ausgedehnt werden kann. Das einfachste Beispiel ist ein 
Kreisring, welches als Scheibe mit einer fast gleich 
grossen Offnung betrachtet werden kann: In diesem Beispiel 
5 sind die Stromlinien 1^ vernachlassigbar (die Stfitzpunkte 
konnten weggelassen werden) und es gibt nur ein 1^ dessen 
Form durch den Kreisring gegeben ist. Falls die Stutzpunkte 
weggelassen werden, ist das Feld B t durch das 
Linienintegral fiber die Gleichung 4 zu bestimmen. 

Ein weiterer Aspekt besteht darin, dass StSreinflfisse von 
unbekannten Objekten abgeschirmt werden, indem man zwischen 
Feldgeneratoreinheit und Objekt eine leitende Platte 
vorsieht, wobei deren Grosse und Form als auch deren 
Position bekannt sind. Damit mussen zwar die 
Feldverzerrungen dieser Platte berucksichtigt werden, doch 
alle anderen elektrisch leitenden Objekte, welche sich auf 
der in Bezug auf die Feldgeneratoreinheit anderen Seite der 
Platte befinden, konnen wegen der Abschirmung 
unberucksichtigt bleiben. 
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Pa tentanspruche : 

1. Verfahren zur Bestimmung der Position eines 
5 sensorelementes (300) , mit Hilfe dessen ein von mindestens 
einer Generatoreinheit (200) ausgesendetes magnet isches 
Wechselfeld (210) gemessen wird, wobei aufgrund eines im 
Sensorelement (300) empfangenen Signals die Position des 
Sensorelementes (300) bestimmt wird, dadurch 
10 gekennzeichnet, 

- dass, vorzugsweise in erster Naherung, Storf elder 
(410) berechnet werden, welche aufgrund von in 
elektrisch leitenden Objekten (400) erzeugten 
Wirbelstromen (420) entstehen, und 

15 - dass die Position, welche ausgehend vom im 

Sensorelement (300) empfangenen Signal bestimmbar ist, 
aufgrund der berechneten Storf elder (410) korrigiert 
wird . 

20 2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, 

- dass die Wirbelstrome (420) im Objekt (400), ausgehend 
vom Wechselfeld (210) , berechnet werden und 

- dass die Storf elder (410) , ausgehend von den 
berechneten Wirbelstrome (420) berechnet werden, 

25 

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass 
zur weiteren Verbesserung der Pbsitionsbestimmung 
mindestens eine weitere Iteration durchgefuhrt wird, indem 
. weitere Wirbelstrome im Objekt (400) , ausgehend von 
30 berechneten Storf eldern (410) , berechnet werden und 

- weitere Storf elder, ausgehend von den weiteren 
Wirbelstromen, berechnet werden. 
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4. Verfahren nach einem der vorangehenden Anspruche, 
dadurch gekennzeichnet, dass Position und Form der Objekte 
(400) bestintmt werden und dass Widerstande nach der 
Gleichung 



im Objekt (400) berechnet werden. 

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass 
15 zur weiteren Verbesserung der Induktivitaten eine weitere 

Iteration durchgefiihrt wird, indem 

- weitere Wirbelstrome im Objekt (400) , ausgehend von 
berechneten Storfeldern (410) , berechnet werden, 

- weiter Induktivitaten, ausgehend von den weit-eren 
20 Wirbelstromen, bestimmt werden. 

6. Verfahren nach Anspruch 4 und 5, dadurch gekennzeichnet, 
dass die Bestiinmung der Objekte (400), die Bestimmung der 
Widerstande und Induktivitaten im Objekt (400) im Sinne 

25 einer Systemkalibrierung im voraus, d.h. vor den das 

Wechselfeld (210) berucksichtigenden Berechnungen, bestimmt 
werden . 

7. Verfahren nach einem der vorangehenden Anspruche, 
30 dadurch gekennzeichnet, dass zur Bestimmung der 




und Induktivitaten nach der Gleichung 
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Wirbelstrome wie folgt vorgegangen wird: 

- die Objekte (400) werden in Segment e (Si) und 
Stutzstellen (Pi) aufgeteilt; 

- in den Stutzpunkten (P A ) wird die Stromdichte (i A ) 
5 bestimmt; 

- aus den Wirbelstromen (1^) werden die Stroradichten 
(i A ) bestimmt. 

8, Anwendung des Verfahrens nach einem der Anspruche 1 bis 
10 7 fur die magnet f eldorientierte Ortung bei „ Cyber Space* - 

Anwendungen . 

9. Vorrichtung zur Durchfuhrung des Verfahrens nach 
Anspruch 1 bis 1, dadurch gekennzeichnet , dass mindestens 
eine Feldgeneratoreinheit (200) , mindestens ein 

15 Sensorelement (300) und eine Verarbeitungs- und 
Steuereinheit (100) vorgesehen sind, wobei die 
Feldgeneratoreinheit (200) und das Sensorelement (300) mit 
der Verarbeitungs- und Steuereinheit (100) verbunden sind. 

10- Vorrichtung nach Anspruch 9, dadurch gekennzeichnet, 
20 dass mindestens ein elektrisch leitendes Objekt zum 

Abschirmen der Feldgeneratoreinheit (200) vorgesehen sind. 

11. Computerprogrammprodukt , das in den internen Speicher 
eines digitalen Computers geladen werden kann und 
Sof twarecodeabschnitte umfasst, mit denen die Schritte 
25 gemass den Anspruchen 1 bis 6 ausgefuhrt werden, wenn das 
Produkt auf einem Computer lauft. 
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